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ABSTRACT 

Diffuse far-ultraviolet stellar emission scattered by dust grains has been observed in a region near 
the Orion Nebula. In addition to the scattered stellar continuum, emission and absorption features 
produced by molecular hydrogen have been identified. In this Letter, we present an analysis of this 
absorption and fluorescent emission from molecular hydrogen in Orion. We model the spectra obtained 
with the Far Ultraviolet Spectroscopic Explorer using optical depth templates and a fluorescent 
emission code. These results are surprising because previous studies have found little ultraviolet 
absorption from H2 in this region, and the emission is coming from a seemingly empty part of the 
nebula. We find that the emission fills in the observed absorption lines where the two overlap. These 
data support the claim that fluorescent excitation by ultraviolet photons is the primary mechanism 
producing the near-infrared emission spectrum observed in the outer regions of the Orion Nebula. 
Subject headings: ISM:molecules — ISM:individual (M 42) — reflection nebulae — ultraviolet:ISM 
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1. INTRODUCTION 

The Orion Nebula (M42) is among the most well- 
studied objects in the sky. Its proximity to the Sun com- 
bined with a bright Photo-Dissociation Region (PDR) 
and a giant molecular cloud make it an ideal region 
to study a range of astrophysical processes. Molecular 
hydrogen (H2) should account for the majority of the 
molecular mass in this region. Observations of emis- 
sion and absorption features of this molecule provide 
useful diagnostics of the physical conditions in PDRs 
and molecular clouds iBurton et al.lll989HLuhm an et al.l 
[ 19941 iRachford et alJ l2002t iKristensen et all 120031 
iFrance et aLl l2005). Previous observations of Ho emis- 
sion in Orion have been limited to the weak quadrupole 
rovibrational transitions arising in the near and mid- 
infrared (IR; see references above as well as Habart et al. 
2004). In this Letter, we present an analysis of the first 
far-ultraviolet fluorescent emission spectra of the elec- 
tronic transitions of H2 in Orion, obtained by the Far 
Ultraviolet Spectroscopic Ex plorer (FUSE ) . Th ese 
spectra have been presented bv lMurthv et alJ l)2005l) . in 
work focusing on the dust scattering properties of the 
region. The spectra show a dust scattered continuum, 
onto which absorption and emission features of H2 are 
identified. 

Molecular hydrogen can be excited into upper rovi- 
brational levels via fluorescent excitation by ultraviolet 
(UV) photons, non-thermal electrons, and collisionally 
(thermally), by shocks. In a gas-rich region that contains 
outflows from young stars and an intense ultraviolet ra- 
diation field produced by O and B stars, it seems likely 
that multiple processes are at work. Near-IR observa- 
tions of H2 in Orion have been used to distinguish which 
of these excitation mechanisms are consistent with the 
obser ved emission ijLuhman et alll994t IKristensen et alJ 
l2003|) . These FUSE observations do not rule out the 
possibility of a contribution from shocks, however they 
provide solid evidence for UV fluorescence in Orion. 

In §2, we present the FUSE data obtained dur- 
ing a serendipitous pointing equilibration observation 



(see iMurthv et ail l|2005[) for a comprehensiv e descrip- 
tion of the observations). IMurthv et alJ l)2005l) find that 
the dust scattered light is dominated by the intense radi- 
ation field of 6 1 Ori C (HD 37022; 06V Maiz-Apellaniz et 
al. 2004), despite the small angular separation of the ob- 
served field from HD 36981 (B5V) . Given its intense far- 
UV output, 8 1 Ori C is likely to be the dominant source 
of excitation for the H2 emission observed in Orion. A 
model of the H2 fluoresence, using the unattenuated con- 
tinuum of 9 1 Ori C as the excitation source, is described 
in §3; that model is compared to the observations in 
§4. These observations provide evidence that UV excita- 
tion gives rise to the observed near-IR emission lines in 
th e Orion Nebula region , in agreement with the findings 
of lLuhman et all l|1994fl . 

2. FUSE OBSERVATIONS AND ANALYSIS 

The diffuse region of the Orion Nebula located at 
RA = 05 ,l 35 m 12.43 s , 5 = -05°lT35.3", J2000, was ob- 
served by FUSE for 16.7 ks in time-tagged (TTAG) 
mode on 2001 November 26 as part of the S405 chan- 
nel rea lign ment program describe d bv lMurthv fc Sahnowl 
l)2004|) and IMurthv et alJ l|2005|) . Nebular spectra were 
acquired across the FUSE bandpass (905 - 1187 A) in 
the (30" x 30") LWRS aperture (Figure 1). Descrip- 
tions of the FU SE instrument and on-orbit performance 
characteristics can be found in Moos et al. (2000) and 
Sahnow et al. (2000). The spectra were obtained from 
the Multi-Mission Archive at STScI (MAST), and were 
analyzed with the original CalFUSE calibration (version 
1.8.7). The calibrated data files for each orbit were then 
coadded using IDL software and the individual channels 
were combined using a cross-correlation algorithm. The 
exact pointing was determined by comparing Fine Er- 
ror Sensor (FES) images with images generated by sky 
simulator software. 

We present the S4054601 spectra at wavelengths longer 
than Ly-(3 (Figure 2), where the bulk of the far-UV emis- 
sion from H2 resides. In addition to the numerous ab- 
sorption lines coming from molecular hydrogen discussed 
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Fig. 1.— 2MASS Atlas image of the Orion Nebula, M42. 1 Ori 
is at the origin and the FU SE LWRS aperture overlay is shown at 
the upper right near HD 36981. The pointing was determined from 
an analysis of the FES image of the target field. Atlas Image mosaic 
obtained as part of the Two Micron All Sky Survey, a joint project 
of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National 
Science Foundation. 



below, we see the scattered photo-spheric lines from 
C III A977/1176, P V A1118/1128, and Si IV A1122/1128, 
as well as interstellar absorption from C II A1036 and 
N I A1135. These features are identified in Figure 2. The 
fact that the P V All 18 line was seen in the spectrum is 
probably the most conclusive argument for 9 1 Ori C (as 
opposed to HD 36981) illuminating this region of Orion. 
Atomic hydrogen lines a re highly contaminat ed by the 
interstellar medium and IPellerin et alJ l|2002D find that 
P V is not present in stars of spectral type later than 
B2V. Finally, we find an unidentified absorption feature 
at 1031.4 A that we tentatively assign to blue-shifted 
O VI. Blue-shifted O VI absorption is seen i n many stars 
of of this spectral type ( Pcllerin et al. 2002) , however we 
cannot conclusively identify it as the A1038 component 
of the doublet falls in a region dominated by H 2 and C II 
absorption. 

3. MODELING MOLECULAR HYDROGEN IN ORION 

The ultraviolet absorption and emission features of H2 
observed in Orion are seen im posed on the dust scat- 
tered continuum of 6 1 Ori C <|Murthv et all 12003) . A 
model of the nebular dust-scattering continuum is be- 
yond the scope of this work, but we set the baseline for 
the molecular hydrogen by fitting the continuum with a 
double power-law of the form 

s{\) = (A)* _ C (±-^r (i) 

Ao Ao 

scaled to the data. This simple emperical expression with 
A of 1187 A, yields a = 0.1, = 2.0, and C = 5.0. 
An H 2 absorption spectrum wa s calculated using the 
H 2 ools optical depth templates l)McCandlissll2003f) . As 
the Doppler 6-value is uncertain, we have adopted 2 km 
s _1 for this diffuse region, consistent with b- values ob- 
served along t he line of sight to other stars in this re- 
gion of Orion ijSpitzer et alJ 119741: Uenkins fc PeimbertJ 



Il997[) . With that assumption, we used the templates 
to fit the (4 - 0) band (A w 1050 A) and find a to- 
tal N(H 2 ) = 4.7 x 10 18 cm" 2 . Absorption out of the 
ground electro- vibrational state for the first six rotational 
states was considered, but most of the absorption comes 
from j = and 1. Because the absoption lines are par- 
tially filled in by the fluorescent emission (see below), the 
true AT (if 2) for this line of sight is most likely greater 
than 4.7 x 10 18 cm~ 2 . This level of H 2 absorption is 
surprising, as studies of stars in this region find a very 
low molecular fraction in the intervening gas (Abel et al. 
2004 and references therein). 

We fit the emission spectrum of H2 seen in these 
FUSE da ta using the hydroge n fluorescence model de- 
scribed in iFrance et alJ l)2005[) . Using parameters for 
the Orion region found in the literature for AT (if 2) and 
temperature, this fluorescent emission model assumes a 
ground electronic state population, then uses photoex- 
citation cross sections and an incident radiation field 
to calculate the rovibrational levels of the upper elec- 
tronic states (predominantly f? 1 ^ and C 1 II U ). The 
molecules then return to the ground electronic state fol- 
lowing the appropriate selection rules and branching ra- 
tios, producing the observed ultraviolet emission lines 
and leaving the molecules in excited rovibrational lev- 
els. iKristensen et alJ £2003) imaged regions of the Orion 
Nebula in the near-IR rovibrational lines of H 2 to de- 
termine the excitation temperature. They find a clear 
delineation between hot and cold zones, requiring more 
than one temperature component to produce a reason- 
able fit to the observations. We ther efore adopted the ro- 
tational temperature of 390 K from lHabart et all l)2004f) 
as the cool (thermal) component and a vibrational tem- 
perature of 2500 K as represe ntative of the hot (s hocked 
or non-thermal) compone nt llFran ce et al] [2 05). An 
N(H 2 ) = 1 x 10 21 cm" 2 (|Habart et al.ll2004D was used 
with a b- value of 2 km s" 1 . Each of these parameters 
affects the resultant model spectrum differently. The 
model spectra change very little with b for values of a 
few km s _1 . Varying the column density changes the 
total emitted power, higher column density giving more 
output photons, but self-absorption begins to surpress 
discrete lines at columns greater than ~ 5 x 10 21 cm~ 2 . 
Temperature controls both the shape and scale of the 
model spectra, with more levels of the ground electronic 
state are populated at higher temperatures. 

Radiation from 8 1 O ri C dominates the ambient field 
in this part of Orion l)Murthv et al.ll2005l) . It is most 
likely responsible for pumping the observed fluorescence, 
and should thus be used as the input field for modeling 
the H2 emission. This is complicated by two factors: 1) a 
well flux-calibrated far-UV spectrum of 8 1 Ori C does not 
exist and 2) the majority of the absorption seen on the 
HD 37022 line of sight is thought to be located in inter - 
vening clouds <)Bautista et al. Ill 9951 iMurthv et al.ll2 005L 
A spectrum of 9 1 Ori C as observed from Earth will not 
be representative of that seen by the molecules in Orion. 
In an attempt to overcome these concerns, we chose a 
lightly reddened SMC star of similar spectral type (AV 
243, 06V, E B - v = 0.09) from the FUSE Magellenic 
Cloud Atlas of Danforth et al (2002; see also Walborn et 
al. 2002). A simple correction for the differences between 
Milky Way and SMC dust was made using the com- 
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Fig. 2. — The FUSE spectrum described in the text overplotted with the H2 emission model described in §3. The model spectrum has 
been convolved to a resolution of 0.5 A. H2 lines are the dominant features seen in both absorption and emission between Ly-/3 and the 
end of the FUSE bandpass. C II and N I are of interstellar origin. The longest wavelength H2 absorption strengths are underpredicted 
due to an offsetting contribution from the emission lines, as described in § 4. 



puted extinction curves of lWeingartner fc Drainei §2001) , 
assuming the reddening of 9 1 Ori C at the FUSE point- 
ing is similar to what is observed along the HD 36981 
sightline (Eb - v — 0.04). We assumed that the red- 
dening on the AV 243 line of sight was due to equal 
contributions from Milky Way [Ry = 5.5) and SMC 
dust. An IUE spectrum of 9 1 Ori C (swp01394, also 
obtained from MAST; Bohlin & Savage 1981) was used 
to cross-calibrate the surrogate FUSE spectrum in the 
overlap region between 1150 and 1187 A. Primary uncer- 
tainties in this approach are a lack of knowledge about 
the attenuation between the star and the nebula, the 
actual contribution from SMC and Milky Way dust, 
and p ossible calibration problem s with the IUE spec- 
trum (Massa fc Fitzpatrickl 12000(1 . Using this incident 
radiation field and the parameters described above, we 
constructed the fluorescence spectrum (convolved with a 
0.5 A box kernel) for comparison with the data. 

4. DISCUSSION 

The nebular spectra and the model described above are 
shown in Figure 2. The model finds satisfactory qualita- 
tive agreement with the relative line strengths seen in the 
FUSE spectra between Ly-/3 and the end of the band- 
pass. We find effects similar to those seen in far-UV 
spectra of the emission/reflection nebula IC 63, namely 



that the observed H2 emission lines are seen to be broader 
than what is expected from the i nstrumental profile for 
these filled aperture observations ((France ct al. 200l|) • A 
study of the possible mechanisms that would produce 
this broadening is underway. The two major discrep- 
ancies that we find between the data and the relative 
strengths predicted by the model are in the region be- 
tween 1040 - 1050 A, and the band of lines centered on 
1161 A. The 1040 -1050 region may be contaminated by 
broad wind emission from O VI, as 9 1 Ori C a lso shows 
an "excess" in this range ijMurthv et alJl2005|) . but the 
exact cause is uncertain. The 1161 band seems to be 
sitting on a "shelf" that is not predicted by the model. 
One hypothesis that we explored was the possibility of 
hot H 2 (T > 2500 K) being fluoresced by nebular Ly- 
a, as i s seen in the far-UV sp ectra of pre-main-sequence 
stars (fWilkinso n el; ftl,j|2QQ2|) and accreting binary sys- 
tems ( Wood & Karovska 2004) . In order to test this the- 
ory, we added a linear continuum that extended the exci- 
tation spectrum to 1300 A and a Gaussian emission line 
at Ly-a. This additional excitation out of excited states 
of H2 does not resolve our problem with the 1161 band. 
We find that as the 1161 band begins to fill in that nu- 
merous other lines are predicted that we do not observe 
in the FU SE spectra. The most likely conclusion is that 
our model is lacking some physical mechanism that gives 
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rise to these features. 

The model does not agree with the observed absolute 
flux. We find that in order to correctly predict the ob- 
served absolute flux, our model output needs to be scaled 
up by a factor of two. We attribute this to three possible 
causes: 1) an uncertainty in the degree of H2 clumping 
in this region, 2) an uncertainty of the true strength of 
the exciting radiation field, and 3) the (unknown) mech- 
anism that causes discrepancies at 1045 and 1161 A. 
Dust and molecular hydrogen are known to be found in 
clumpy and/or filamentary structures in r egions where 
hot stars interact with their surrounding s (O'Dcll1 l2000t 
iKristensen et al.l 120031 iFrance et all 12004ft . however we 
find this explanation unlikely. Creating similar models 
at higher column densities shows that while the overall 
emitted power increases with column, the majority of this 
extra emission is located outside of the FUSE bandpass 
(A > 1200 A). The second possibility is that our input 
excitation field does not accurately represent the local 
absolute radiation field. We feel confident that the spec- 
tral characteristics of the input spectrum are good based 
on the relative agreement with the H2 emission spectrum, 
but are less confident in the absolute scale. Changes in 
our reddening assumpti ons, absolute calibration o f the 
IUE spectrum we used (Massa fc Fitzpatricldl200C)|) . and 
uncertainties in the nebular distance/geometry can easily 
explain the scaling required. Finally, further modeling 
efforts are underway to resolve the problems described 
above. 

The molecular hydrogen absorption lines present 
in the spectra also show an interesting behavior. 
One woul d expect, based o n the decreasing oscillator 
strengths ( Mc Candliss! 120031) . that the lines would be 
shallower during the progression to smaller vibrational 
transitions (4 - to - 0, shown in Figure 2). We find 
that the absorption profiles are of almost identical depth 
in each of the bands. We attribute this to a filling ef- 
fect by the emission from H2 where the emitting and 



absorbing transitions overlap. There are a greater num- 
ber of emitting transitions overlapping with the absorp- 
tion lines at shorter wavelengths, so fitting the absorption 
lines at the (4 - 0) band leads to an under prediction of 
the strength of the lines in the longer wavelength bands 
of molecular hydrogen. This finding leads us to stress 
that meaningful column densities cannot be determined 
without correcting for the emission lines and properly ad- 
dressing the dust scattering to set the continuum level. 

There has been uncertainty in the degree to which 
shocks or ultraviolet photo-excitation leads to the 
observed n ear-IR lines o f mo le cular hydrogen seen 
in Orion llBurton et alJ 119891 iLuhman etahl Il994t 
IKristensen et alJ20 031. These FUSE data rule out mod- 
els that consider only shocks as the mechanism produc- 
ing the near-IR emission. Shocks cannot populate the 
excited electronic states of mol ecular hydrogen that giv e 
rise to the emission lines seen (IShull fc Beckwitbl ll982). 
Hydrogen molecules will be dissociated by collisions at 
temperatures greater than roughly 4000 K, thus they 
would be destroyed before reaching temperatures suf- 
ficient to begin populating the first excited electronic 
state, -B 1 ^ (E « 11 eV). These data cannot exclude 
shocks fr om making a contribut ion to the near-IR H2 
emission l|Kristensen et all 120031) . but ultraviolet exci- 
tation clearly contributes to the rovibrationally excited 
population of hydrogen molecules in the Orion Nebula. 
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